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RADIATION  FROM  INTENSE  ELECTRON  BEAMS  ASSOCIATED  WITH 

THE  CERENKOV  MECHANISM 
John  R.  Neighbours  and  Fred  R.  Buskirk 

Abstract 

This  report  summarizes  the  Cerenkov  radiation  to  be  expected 
from  an  intense  electron  beam  for  both  the  coherent  and  incoherent 
modes  of  radiation.  The  effects  of  a  finite  length  beam  are  dis- 
cussed with  regard  to  oeam  diffraction,  radiated  power,  and  bunch 
snape.  The  appendices  contain  a  discussion  of  the  beam  radiation 
pattern . 


RADIATION  FROM  INTENSE  ELECTRON  BEAMS  ASSOCIATED  WITH 

THE  CERENKOV  MECHANISM 
1.   INTRODUCTION 

Normally  a  charge  or  system  of  charges  moving  at  constant 
velocity  does  not  radiate  electromagnetic  energy,  but  if  the 
velocity  of  the  charge  is  greater  than  the  velocity  of  light  in  a 
medium,  radiation  occurs.   It  is  broad  band;  the  intensity 
increases  with  frequency;  and  the  propagation  is  confined  to  a 
cone,  much  like  a  shock  wave  in  air.   Equations  {?.)    and  ( 3 )  in  a 
later  section  describe  these  results.   In  previous  work  l~ 4 ,  we 
described  theory  and  experiments  showing  that  this  mechanism 
produces  microwaves  from  bunched  electron  beams,  in  addition  to 
the  usual  optical  radiation,  as  studied  by  Cerenkov  and  many 
others . 

In  this  report  we  consider  aspects  of  coherent  Cerenkov 
radiation,  with  emphasis  on  radio  and  microwave  frequencies.   As 
will  be  shown,  power  levels  may  become  very  high,  in  the  megawatt 
range,  so  that  serious  consideration  should  be  given  to  this 
mechanism.   The  scope  of  this  report  is  to  outline  the  systematics 
of  a  theory  in  which  the  dielectric  (v/hich  lowers  the  velocity  of 
light)  is  considered  to  be  linear  and  non-dispersive.   The  aspects 
considered  which  are  not  treated  in  the  extensive  literature  are 
calculations  for  electrons  occurring  in  finite  size  bunches,  both 
single  and  repeated.   Future  work  should  aim  at  the  effect  of 
ionization  in  the  dielectric  medium,  betatron  oscillations  of  the 
beam  electrons,  and  other  mechanisms  which  could  seriously  modify 
or  degrade  the  radiation. 


2.   SUMMARY  OF  PREVIOUS  WORK 

In  the  earlier  work1-4,  it  was  assumed  that  electrons  were 
emitted  from  an  accelerator  in  bunches,  described  by  a  charge 
density  p(r),  and  all  had  the  same  velocity  v  such  that  the  bunch 
did  not  change  shape  in  time.   The  current  density  is  then  given 
by 

J(r)  =  p  v  (1) 

A.   Intensity 

Under  these  conditions,  the  Cerenkov  radiation  in  a  medium 
with  the  velocity  of  radiation  c  <  v  is 


dE       y_ 

dx       4  it       ""  ^c 


ojdw  sin20   q2  F2(k)  (2) 


dE        jj 


<j)'ji 


sin23  a2  F2(k)  (3) 


dx        4  TT      0        c 

Eq.  (2)  refers  to  a  single  bunch,  the  radiation  has  a 
continuous  frequency  spectrum,  and  the  formula  gives  the  energy 
radiated  at  angular  frequency  w  in  a  range  dm,  per  unit  length. 
Eq .  (3)  is  the  energy-per-unit  path  per  bunch  for  the  case  of 
periodic  bunches,  emitted  at  an  angular  frequency  w0 •   Here  w  is 
a  harmonic  of  oo0,  in  contrast  to  the  continuous  frequency 
spectrum  of  Eq.  (2) .   In  both  cases,  q  is  the  total  charge  of  one 

bunch,  and  F(k)  is  a  form  factor  for  the  charge  distribution 

,  -*■    •*■ 
F(k)  =  q"1///  p(r)  elk'r  dV  (4) 

0C  is  the  Cerenkov  angle, 


q  is  the  charge  per  bunch,  and  k  is  propagation  vector  for  the 
emitted  radiation,  such  that  k  =  co/c,  and  the  direction  is  alonq 
the  Cerenkov  cone. 

B.  Polarization 

The  electron  field  is  in  the  plane  of  the  observer  and  the 
electron  path,  as  it  is  for  usual  optical  Cerenkov  radiation. 

C.  Coherence 

If  the  charqe  is  a  point  charge,  F  =  1  and  in  fact,  for  a 
finite  bunch,  we  will  have  F  ~  1  for  \l    <<    1,  where  k  =  2  it/  X . 
This  corresponds  to  X  >  %,    where  I    is  the  size  parameter  of  the 
bunch.   For  these  long  wave  lengths,  either  formula,  Eq .  (2)  or 
(3),  shows  that  the  intensity  is  found  by  setting  q  =  ne  where  n 
is  the  number  of  electrons  in  the  bunch.   Thus  the  intensity  is 
proportional  to  n2e2.   As  higher  frequencies  are  examined,  the 
factor  F  becomes  small,  and  finally,  for  very  high  frequencies 
(optical  range),  we  would  expect  incoherent  radiation,  in  which  q2 
is  replaced  by  ne2. 

P.   Diffraction  Effect 

An  important  result  is  that  if  the  region  of  emission  of 
radiation  is  finite,  of  length  L,  the  emitted  radiation  is  no 
longer  confined  to  the  Cerenkov  cone  0  =  9C  but  is  spread  about 
that  angle  and  the  spreading  depends  on  the  length  of  the  path  in 
the  gas  and  the  wave  length  of  the  radiation.   It  is  more  closely 
related  to  the  diffraction  pattern  of  an  end  fire  antenna,  rather 
than  the  broadside  array;   the  latter  is  related  to  the  single 
slit  of  optics.   This  effect  is  independent  of  bunch  size  and 
repetition  rate.   Associated  with  diffraction  is  an  enhancement, 
qualitatively  described  below. 


E.  Enhancement  from  Diffraction 

In  most  diffraction  or  interference  situations  in  optics, 
the  energy  is  usually  redistributed  in  angle  if  parameters  are 
changed.   For  example,  if  the  width  of  a  single  slit  is  changed 
but  the  total  power  passing  through  the  slit  remains  constant,  the 
total  power  integrated  over  the  diffraction  pattern  on  a  screen 
remains  constant.   In  the  Cerenkov  case,  the  development  of  Ref.  3 
showed  that  for  small  Cerenkov  anqles,  the  intensity  at  0C 
remains  constant  as  L  decreases,  but  the  total  power  radiated  at 
angles   0  >  9C  increases  dramatically.   The  total  power  radiated 
increases  by  factors  of  up  to  50  in  our  experimental  situation. 

F.  Form  Factor  Effect 

The  basic  Cerenkov  radiation  strength  tends  to  increase 
with  frequency,  for  either  coherent  or  incoherent  emission. 
Starting  at  lov;  freauencies,  the  coherent  emission  increases  with 
co  until  F  becomes  small,  at  oj  =  u^.   At  this  frequency,  the  wave 
length  is  short  enough  so  that  radiation  emitted  from  electrons  at 
the  front  and  rear  of  the  bunch  begin  to  interfere  destructively. 
This  leads  to  destructive  interference  and  lowered  emissions  of 
power  above  co^.   Experiments  in  principle  could  demonstrate  this 
effect,  but  at  this  time  we  have  been  unable  to  make  measurements 
in  the  high  frequency  microwave  range  (  A  ~  2  mm)  which  would 
quantitatively  demonstrate  the  predictions. 

As  a  further  remark,  it  should  be  possible  with  precise 
measurements  to  determine  size  parameters  of  the  electron  bunches 
at  a  distance,  or,  in  other  words,  to  do  beam  diagnostics  by  this 


method.   These  possibilities  will  be  explored  in  another  report. 
3.   Approximate  Results  for  Radiated  Power 

The  results  mentioned  previously  may  be  related  to  more 
conventional  parameters  of  pulse  length  and  peak  current.   If  we 
assume  a  pulse  in  the  form  of  a  uniform  cylinder,  J  =  pv  and  the 
current  I  =  JA,  q  =  pAl    =  pAvT  where  I   =  the  length  of  the  bunch 
and  T  =  the  bunch  length  measured  as  a  time.   Thus  we  find 

q  =  IT  (6) 

For    a    Gaussian    pulse,    the    charge   density    is 

-(x2   +   y2)/a2      -b2z2/b2  (7) 

p    =    p0    e  -r  e 

It    may   be    shown    that    the    peak    current    in    the   middle    of    the 

pulse    (at    z    =    0)     is 


I       =    -^—    q  (8) 

o  Tib 

and  that  the  normalization  constant  P   is  related  to  the  total 


charge  of  the  bunch  by 


3/2    2 
p   it     a  b  (9) 


Then  we  may  let  the  effective  length  of  the  bunch  be  I   =   bir 

and  use  q  =  L  T  =  L  £/v,  just  as  for  the  uniform  cylinder. 

Now  consider  the  sum  over  frequencies  in  Eq.  (3).   Assuming 

that  F(k)  is  constant  up  to  some  k     which  corresponds  to 

max  r 

a)   =  Nojn,  we  find 
m      u 


5>u)Q    =    £nu)02 


=    N(N+1)     wn2/2 


N2wn2/2 


0  ^liUJo 

=    -m2/2 
Thus 

£u>u>0    =    0^2/2  (10) 

For  a  single  bunch,  the  sum  over  the   continuous  frequency 
distribution  (2)  gives 

f     mU)d03     =      U3m2/2  (11) 

Thus  the  results  of  Eqs .  (10)  and  (11)  above  may  be  used  for 
either  single  or  repeated  pulses  to  calculate  the  total  power 
radiated  at  all  frequencies. 

Now  we  may  use  the  results  of  Eq.  (6)- (9)  to  re-express  the 
results  of  Eq.  (?)  or  (3).   The  latter  may  be  multiplied  by  v  to 
convert  to  power,  then  summed  or  integrated  over  oj  for  the  total 
power  emitted,  which  is 

2    2    2 

P  =  ^0-  v  sin20   In  T  f  (12) 

2  c   u     m 

This  holds  for  either  periodic  or  a  single  bunch.   In  either 
case  it  is  the  rate  at  which  electromagnetic  radiant  energy  is 
emitted  while  the  bunch  is  in  flight,  per  bunch  (in  the  case  of 
periodic  bunches). 

Noting  that  fm  =  l/T  holds,  which  is  equivalent  to  the 
uncertainty  relation  in  quantum  theory  or  the  bandwidth-pulse 
lenqth  relation  in  communications, 


P  =  JDLq-   v  sin20   I02  (13) 


Thus  the  power  emitted  for  coherent  radiation  is  proportional 
to  the  square  of  the  beam  current.   This  makes  explicit  one  of  the 
main  points  of  this  report:  that  the  emitted  power  increases 
rapidly  with  the  beam  current.   Thus,  although  weak  beams  emit 
little  radiation,  the  radiation  becomes  important  for  large 
currents.   When  we  note  that  the  radiation  emitted  is  proportional 
to  Iq2  but  the  energy  in  the  beam  is  proportional  to  IQ,  so  that 
relative  efficiency  as  a  radiator  increases  in  proportion  to  IQ. 
This  may  be  made  quantitative  by  multiplying  (13)  by  L/v  to  aet 
the  total  energy  emitted  by  the  bunch  during  its  flight  through  a 
path  L,  where  L  is  the  range.   The  energy  put  into  the  bunch  is 
IqT,  so  that  the  efficiency  e,    defined  as  the  ratio  of  the  energy 
radiated  (Cerenkov)  to  the  total  energy  of  the  beam  becomes 

,  =  If-   Sin\  -i*  JL.  (14) 

where  VQ  is  the  energy  of  the  beam  expressed  in  volts  and  L  is  the 
range  of  the  beam  in  air.   Strictly  speakinq,  L  is  the  distance  at 
which  Cerenkov  radiation  ceases.   The  precise  value  of  L  would 
depend  on  the  relative  importance  of  bremsstrahlung,  ionization, 
and  instabilities  which  could  disperse  the  electron  bunch,  and 
these  aspects  should  be  considered  in  greater  detail . 

The  overall  systematics  of  the  radiation  associated  with  the 
Cerenkov  mechanism  are  expressed  by  Eq.  (12),  (13),  and  (14).   The 
following  features  should  be  noted: 


A.  The  power  radiated  per  bunch  (eq.  12,  13)  is  proportional 
to  the  square  of  the  peak  beam  current,  ^q2'  arv^  ^s  independent  of 
the  pulse  length  T.   This  depends  on  the  approximation  that 

fmT  =  1,  where  T  is  the  pulse  length  and  fm  is  the  maximum 
frequency  for  which  the  fourier  transform  of  the  pulse  has 
significant  strength.   This  expression  for  the  power  radiated  is 
the  energy-per-unit  time  while  the  bunch  is  in  flight,  from  the 
accelerator  to  the  point  where  it  ceases  to  radiate  effectively. 
The  latter  point  in  the  simplest  consideration  is  the  point  at 
which  the  bunch  has  lost  enough  energy  so  that  it  is  below  the 
threshold  for  Cerenkov  radiation,  which  is  20  MeV  for  electrons  in 
normal  air . 

B.  As  T  is  decreased,  assuming  IQ  is  constant,  the  power 
radiated  remains  constant.   But  the  frequency  at  which  most  of  the 
radiation  occurs  is  increased  as  T  is  decreased. 

C .  The  efficiency  (Eq.  14)  depends  on  I0/T  showing  that  for 
high  currents  and  short  bunches,  the  Cerenkov  mechanism  could  not 
only  lead  to  sizable  amounts  of  radiation  but  also  dominate  the 
other  modes  of  energy  loss. 

D.  As  a  guide  to  what  may  be  expected  in  the  case  of  plasma 
shielding,  we  might  expect  that  some  length  I'    of  the  nose  of  the 
beam  bunch  might  be  unshielded  by  the  plasma.   From  the  above 
considerations,  the  length  V    would  not  change  the  total  power 
radiated,  but  would  change  dramatically  the  freauency  at  which  the 
radiation  occurs.   This  is  because  if  only  the  nose  region  of 


of  length  V    radiated,  the  total  effective  radiating  charge  would 
be  proportional  to  % ' .   But  the  frequency  of  the  maximum  radiation 
would  increase  as  1/ i' ,    so  the  power  radiated  would  remain 
constant  as  V    is  varied. 
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APPENDIX  A  -  EXTREME  VALUES  OF  THE  RADIATION  FUNCTION 

The  radiated  power  per  unit  solid  angle  W(v,n)  can  be  written  as  the 
product  of  two  factors: 

W  =  Q  R2  (Al) 

where  Q  is  a  function  having  the  dimensions  of  power  (watts)  and  R  is  a 
dimensionless  radiation  function. 

Q  =  q2  (A2) 

8tt2 

R  =  (kL)  SinG  F(k)  I(u)  (A3) 

Here,  u  is  the  permeability  of  the  medium,  c  is  the  velocity  of  light  in 

the  medium,  and  q  is  the  charge  of  an  individual  electron  bunch.   If  an  RF 

Linac  is  the  source  of  the  electrons,  the  frequency  vQ  is  that  of  the  Linac ; 

however,  in  all  cases,  vQ  corresponds  to  the  spatial  periodicity,  I  =   c/v0, 

of  the  electron  bunches. 

In  (3),  k  is  the  wave  vector  of  the  emitted  radiation,  L  is  the  beam 

interaction  length,  and  F(k)  is  the  form  factor  defined  such  that 

p0'(k)  =  q  F(k)  (A4) 

and  £Q'(k)  is  the  Fourier  transform  of  the  charge  distribution  of  an  individual 

charge  bunch.   The  diffraction  function 

I  =  Sln  U  (A5) 

u 

is    a    function    of    the    diffraction   variable 

u  =   kL    (Cos9     -    CosOl  (A6) 

"2 

and  consequently,  R  is  a  somewhat  comDlicated  function  of  0  and  kL . 
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MINIMUM  VALUES 

By  inspection,  it  is  obvious  that  R2  and  thus  W  will  have  minima   when 

l(u)  =  0.   These  true  zeroes  of  the  function  occur  when  u  is  equal  to  integral 

values  of  tt;  i.e.,  u  =  mir.   Substituting  this  condition  in  (A6)  leads  to  the 

angular  values  for  the  zeroes  in  the  radiated  power 

Cos  9  =  Cos  3  -  m  ~  (A7) 

c       kL 

which    can   be    written    in    terms    of    the   wavelength    (in    the   medium)    of    the 
radiation    as 

CosO  -   CosO  -  m    [£)  (A8) 

This    relation    is    exact    and    gives    the    zeroes    of   R   for    all    allowed   values    of 
m   for   which    CosO    <    1.       I(u)    is    a  maximum   at   m   =   0    (which    gives   u   =    0),    and    the 
maximum    in    R  tends    to    occur    somewhere    in    that    vicinity.      The    first    zero    is    that 
one   which    is    adjacent    to    the   maximum   and    is    given   by    (A8)    when  m   =   +1.      Values 
of    this    angle    as    a    function    of  —     are    shown    as    a   Log-Log    plot    in    Fig.    1.  The 

A 

relationship  between  9  and  -r  is  approximately  that  of  the  inverse  square  root 
over  the  range  10  <  L/ X  <  100.   Within  this  range,  the  angular  value  (in 
degrees)  of  the  first  minimum  is  given  bv  the  empirical  relation 

01  =   ,82  (A9) 

VlTx 

The  physical  range  of  9  is  from  0  to  tt.   Substitution  of  the  successively 

larger  positive  values  of  m  =  2,  3,  U...    into  (A7)  or  (A8)  gives  a  set  of 

increasing  angles  for  the  functional  zeroes.   If  X/L  <<  1,  many  zeroes  will 
occur  before  the  physical  limit  is  reached. 

In  addition,  if  X/L  <<  1,  there  is  the  possibility  of  having  functional 

zeroes  at  angular  values  less  than  that  of  the  main  peak.  These  occur  with  the 
substitution  of  successive  negative  m  values  into  (A7)  or  (A8),  giving  rise  to 
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subsidiary  radiation  maxima  between  the  beam  line  and  the  main  diffraction 

peak. 

MAXIMUM  VALUES 

To  find  the  maximum  values  of  Wr  -*  one  proceeds  in  the  usual  fashion  by 
taking  the  derivative  and  setting  it  to  zero. 


||=2QR-|=0  (A10) 

3R 
This  equation  is  satisfied  for  either  R  =  0  or  -r^  =  0.   The  first  condition  is 

that  considered  above  for  the  zeroes  of  the  function.   The  remaining  one, 

3R 
-~7r  =  0,  pertains  to  the  maxima.   Thus,  it  is  sufficient  to  consider  maximum 

values  of  R  which  will  also  be  maximum  values  for  W  ~  R  . 

From  (A3),  the  derivative  of  R  is 

||  -  kL  {  sinO  [F  -||  +  I  -||]  +  FI  Cos9  }  =  0  (All) 


and  the  derivative  of  the  diffraction  function  I(u)  can  be  written 

31  _  dl  3u  _  _,  kL  _.  n  /.nv 

"39  "  du-70  "  l      2~   SinG  (A12) 


Substituting  this  value  in  (All)  and  discounting  the  solution  for  kL  =  0  gives 


the  condition  for  maxima 

kL  _. 
—  Sin. 


{  SinO  [!'  F  ~  SinO  +  I  4^j  +  FI  CosO}  =  0  (A13) 


In  order  to  proceed  further,  the  functional  dependence  of  the  form  factor 
F  must  be  known. 
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A  gausian    form    factor   has    been    assumed    in    this    and   other    reports. 

F(k)    =   Exp    (-a2  Cos29)  (A14) 

where  a=  kb/4  corresponding  to  a  one-dimensional  charge  distribution 

Z2 
PO*  (r)  =  q  Exp  ( )  (A15) 

b2 


The  derivative  of  this  form  factor  is 

||  =  a2  Sin  29  F  (A16) 

which  is  small  for  the  parameters  of  interest  in  Linac  experiments;  i.e., 
b  =  0.24  cm,  microwave  frequencies,  and  angles  up  to  approximately  25°. 
Therefore,  a  first  approximation  to  solving  (A13)  can  be  obtained  by  setting 
the  derivative  of  the  form  factor  in  (A13)  to  zero.   The  resulting  aDproximate 
maxima  condition  is  independent  of  F  and  can  be  rearranged  so  that  I  and  its 

derivative  are  on  the  same  side. 

I'(u)  _  2   Cos  9  ..... 

2-5  (A17> 


I(u)    kL  Sinz  9 
The  ratio  of  I'  to  I  can  be  expanded  as  a  power  series  in  u  as 

which  converges  rapidly  for  small  values  of  u.   Substituting  the  first  term  in 
the  expansion  into  (A17)  and  subsequently  eliminating  the  diffraction  variable 
through  its  definition  (A6)  leads  to  an  approximate  equation  for  the  maxima, 
which  is  cubic  in  y  =  Cos9. 

-Y3  -  Cos9   y2  -  [l  +  — ]  y  +  Cos9p  =  0  (A19) 

C         L       n  J  C 

X 

In   this    equation,    the   beam   parameters    enter    through    the    charge  velocity    in 

Cos    9C    =   c/v    and    the    radiation    wavelengths    in   x   =   kL   =    2 it  L/X.  Numerical 

solutions    can   be    carried    out    on    a  hand    calculator   with    results  within 

approximately   0.3°    of    those    obtained    from    large-scale    computer  calculations. 
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Results    for    the   major    radiation    lobe    (the    first   maximum)    are    also    plotted    in 
Fig.    1.      These    also    show   a    straight    line    dependence    on    a   Log-Log    plot    and 

therefore    can   also   be    represented    as    a  power    law 
1 


0     =52      .__ 
which    is    accurate    to    approximately    5% 
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APPENDIX  B  -  EFFECT  OF  BEAM  LENGTH  ON  THE  RADIATION  FUNCTION 

In  previous  reports,  several  calculated  curves  of  a  radiation  function 
D,  _v  which  differs  only  by  constants  from  R2  were  presented.   Fig.  2  is 
reproduced  from  Fig.  5  of  Reference  4  and  relabeled  to  be  consistent  with  the 
notation  of  Appendix  A.   Fig.  2  shows  the  vartiation  of  the  radiation  function 
for  8.55  GHz  radiation  from  the  MPS  Linac  for  beam  interaction  lengths  in  the 
vicinity  of  100  cm.   As  noted  previously  in  this  and  other  reprots,  increasing 
L  tends  to  move  the  main  diffraction  lobe  toward  the  beam  line,  and  Appendix  A 
gives  an  approximate  formula  for  calculating  the  position  of  the  first  lobe. 
Fig.  3  (to  different  scales)  shows  this  effect  for  larger  beam  paths  from  an  S- 
Band  Linac.   As  L  increases,  the  principle  lobe  moves  in  toward  the  Cerenkov 
angle  3C ;  i.e.,  0^    ■*   0C ,  and  the  radiated  power  per  unit  solid  angle  at 
the  maximum  increases  greatly. 
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FIGURE  CAPTIONS 

Fig.  1   Cerenkov  radiation  extremum  angles  in  degrees  vs  L/A.   The  wavelength  X 
of  the  Cerenkov  radiation  is  measured  in  the  medium  in  which  the  path 
length  is  L.   Curve  B  is  for  the  first  minimum,  the  solution  to  (A8) 
with  m  =  +1 .   Curve  A  is  for  the  first  maximum,  the  solution  to  the 
approximate  equation  (A19).   Over  the  range  shown,  A  is  straight  and  B 
is  slightly  curved.   Equations  for  the  curves  are  given  in  the  text. 

Fig.  2   Calculated  third  harmonic  (8.55  GHz)  radiation  patterns  produced  by 

electron  bunches  from  an  S-band  (2.85  GHz)  Linac .   The  path  lengths  in 
air  are  70,  90,  and  150  cm.   Vertical  scale  is  arbitrary  and  the  same 
for  all  three  curves. 

Fig.  3   Calculated  third  harmonic  radiation  patterns  produced  by  electron 

bunches  from  an  S-band  Linac.   The  solid  curve  is  for  a  path  in  air  of 
1000  cm,  and  the  dotted  curve  is  for  a  path  length  of  10,000  cm. 
Vertical  scale  is  arbitrary  and  the  same  for  both  curves  but  much 
reduced  from  that  of  Fi?.  2. 
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